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We present a microwave realization of finite tight-binding graphene-like structures. The structures 
are realized using disks with a high index of refraction. The disks are placed on a metallic surface 
while a second surface is adjusted atop the discs, such that the waves coupling the disks in the air 
are evanescent, leading to the tight-binding behavior. In reflection measurements the Dirac point 
and a linear increase close to the Dirac point is observed, if the measurement is performed inside 
the sample. Resonances due to edge states are found close to the Dirac point if the measurements 
are performed at the zigzag-edge or at the corner in case of a broken benzene ring. 
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Due to its electronic shell structure carbon forms struc- 
tures such as diamonds, fullerenes and nanotubes, which 
all have very specific and fascinating properties, mechan- 
ically as electronically. Another realization is graphene, a 
one-atom-thick allotrope of carbon, which has a structure 
of the honeycomb lattice, i. e. two combined triangular 
lattices. A common feature of all sheet-like carbon struc- 
tures including graphene is sp 2 hybridization leaving one 
p z orbital for the bonding. This is the starting point for 
the tight-binding approximation used in many theoreti- 
cal descriptions of graphene including the disseminating 
works by Wallace^ and Semenoffj 2 - Due to the specific 
symmetry of the lattice the gapless band structure has 
a conical singularity at two k vectors^ called nowadays 
the Dirac points. The name comes from the fact that the 
reduced equation around the Dirac points corresponds 
to the Dirac equation of a massless particle resulting in 
a band structure of the form of two cones touching at 
the tips.— & Recently it was shown that this single layer 
can be realized and investigated experimentally£i£ The 
main ingredient for this relation is the symmetry. Thus 
a realization via a photonic crystal^ was proposed^ In 
addition in two microwave experiments signatures of the 
Dirac point had been found in the transmission in an 
triangular set-up i 10 ' 11 Both experiments correspond to 
open scattering arrangements. Graphene, on the other 
hand, corresponds to a closed system because of the tight- 
binding between nearest neighbors preventing an escape 
to the outer world. We consider the aspect of symme- 
try and the aspect of tight-binding with mainly nearest 
neighbor interaction equally important.! 2 This is realized 
by a microwave experiment using dielectric disks with a 
high index of refraction in a hexagonal lattice. The disks 
are weakly coupled by evanescent modes as will be ex- 
plained later. 

So far microwave experiments have been used to show 
different kinds of effects originally coming from con- 
densed matter physics such as the Hofstadter butterfly jl2 
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FIG. 1. (color online) Spectrum of a single disk for the TE 
antenna (dark, black) and the TM antenna (light, yellow). 
The TEi resonance of the disc, which is used for further in- 
vestigations is marked. The inset shows the TE antenna with 
a single disk on the metallic copper plate. The top plate is 
not shown. 

localization in disordered systems ; 14 ' 15 and transport in 
case of correlated disorder^ 

In the proposed microwave setup it is feasible to in- 
vestigate defects, edges, disorder etc., as the realized 
"graphene" is a macroscopic object. Thus it is possible to 
study problems experimentally which are not accessible 
in real graphene at the moment. Additionally one can 
have a look into the system and study e. g. the density 
of states (DOS) or transport properties as well. In this 
paper we will concentrate on the local DOS (LDOS) and 
its characteristics at edges or corners. 

As a starting point we investigate a single disc, then 
the coupling of two of them and proceed via a benzene 
ring to graphene flakes. We measure the reflection S\\ of 
an antenna situated close to a disk of r<j=4 mm radius, 
hd=5 mm height and an index of refraction rid ~ 6. The 
antenna is bended (see inset of Fig. [T|) to excite both 
transverse magnetic (TM) and transverse electric (TE) 
modes whereas a standard dipole antenna can only ex- 
cite the TM modes. The disk is placed on a metallic 
plate (see inset of Fig. [T]), a second plate is covering the 
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FIG. 2. (color online) Splitting of the eigenfrequencies of two 
disks as a function of the distance caused by evanescent cou- 
pling. The experimental data (crosses) can be described by a 
Bessel function Ko (see Eqs and (J3J)) where ko ~ 1.82 GHz 
and 7 w 0.288 mm" 1 and a small offset Av » 0.022 GHz. 
The inset shows the spectrum and splitting for d=3.3 mm. 

whole setup at height h = 16 mm. The same disks have 
been used in Ref to find localization in disordered mi- 
crowave structures. In contrast to previous experiments 
the top plate is not touching the disks but it is 11 mm 
above the disc. The distance is chosen such that the reso- 
nance is still sharp but the eigenfrequency is not sensitive 
on small deviations of the distance to the top plate. In 
Fig.[T]the measured spectrum, more precisely 1 — |Sii| 2 , 
of a single disk is shown. The dark black curve presents 
the spectrum measured by the antenna shown in the in- 
set coupling both to the TE and TM modes whereas the 
light curve was measured with a pure dipole antenna only 
coupling to the TM mode. The TEoi -resonance which 
is used solely for the further investigations is marked. 
For all measurements we took care that the resonances 
below and above the investigated frequency range were 
still separated from the frequency range of interest. The 
eigenfrequency of the TE resonance decreased with in- 
creasing height of the plate, demonstrating that it is a 
TEcu-resonance of the disc. 

Now we investigate the coupling of two discs. Be- 
tween the disks the TEi mode can be described approx- 
imately by only a z component of the magnetic field 
B = (0, 0, B z ) and perpendicular components of the elec- 
tric field. Let us assume for the moment that both the 
bottom and top plates touch the disks from above and 
below. Then the z dependence can be separated, yielding 
for the lowest TE-mode close to an individual disk 



B z (x,y,z,k) = B 



sin (jk z ) J oi k ± r ) r <r d 
r > r d 



asm [fz) K Q (^r) 



(1) 

where k± = \J k 2 — (^) 2 \ 7 = \J (f) 2 — k 2 , and r is the 
distance from the center of the disc. Jo and Ko are Bessel 
functions and a is a constant to be determined from the 
continuity equations at the surface. At the resonance fre- 
quency of the disk the wave number for r > r d is thus 
purely imaginary leading to an evanescent coupling be- 
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FIG. 3. (color online) Spectrum of a single benzene cell. The 
dotted line corresponds to a spectrum of a perturbed benzene 
cell, where the lifting of the degeneracy can be seen. 

tween the discs. Since in reality there is a gap between 
the top plate and the discs, this can be correct only ap- 
proximately. But still we use Eq. ([I} with an effective 7 
to be determined from the experiment. A direct coupling 
to the TMo and TMi modes is strongly suppressed due 
to the continuity conditions of the fields. From Eq. (pj 
we obtain for the coupling between two disks 



At(d) = Kq 



#0 7 



(2) 



where d is the distance between the disks (d=0 means the 
disks are touching) and kq is a constant. The frequency 
splitting is given by 



Av(d) = ^4K(d) 2 + (Av ) 2 , 



(3) 



where Ai>q is taking into account a small eigenfrequency 
difference of the disks and the differing coupling of the 
antenna to the discs, which actually dominates. The res- 
onance splitting is shown in Fig. [5] exhibiting an agree- 
ment with the predicted behavior. Thus each disk brings 
in just one bound state coupled evanescently to its neigh- 
bors, exactly the situation found in a tight-binding sys- 
tem. 

As a next test of the tight-binding approximation, we 
performed a measurement with three disks in a row as a 
function of the distance and found the coupling constant 
to be consistent with the two disk measurement. The 
next-nearest-neighbor coupling showed up to be about 
6% of the nearest neighbor coupling. As graphene is 
made of coupled benzene cells we investigated the ben- 
zene cell (see inset of Fig. [3]) next. In Fig.[3]the spectrum 
(dark solid line) is shown. Benzene has dihedral symme- 
try D^h with two singlets and two doublets, where the 
singlets have the extremal energies. Correspondingly we 
observe four resonances. The degeneracy of the doublets 
can be lifted by a perturbation, e. g., by moving one disk 
(see dotted line). Again we performed a measurement 
varying the distance and found a good agreement with 
the coupling constant obtained by the two disk measure- 
ment. 



3 



0.8 




n n I , , , , i , , , , i , , , ; i , , , , i , , , , i , , , , i , 

6.50 6.55 6.60 6.65 6.70 6.75 6.80 
1//GHZ 



FIG. 4. (Color online) Spectrum of a graphene flake, where 
the spectrum was averaged over three different flake forms and 
different interior antenna positions. The dashed vertical line 
marks the resonance position of a single disk corresponding 
approximately to the Dirac point. The dotted horizontal line 
corresponds to the background. The dashed lines give a guide 
to the eye to see the linear increase of the DOS. 

Before we come to the final measurements we have to 
establish a connection between our measured quantity, 
the reflection Sn, and the DOS or more precisely the 
LDOS. The starting point is the relation between the 
Green's function and the local density of states C(E) = 
— -ImG. 17 Here we illustrate a short hand argument the 
connection between the reflection and the LDOS. The 
total current through the antenna is given by T = 1 — 
|<Sii| 2 — § ds j, where j = Im^'*V v E' is the current density. 
The wave function of the field excited by a point-like 
antenna at R is proportional to the Green's function, 

j ds j = i I ds (<fw* - V*) 

= -\B(R;E)\ 2 lmG(Il;E). (4) 

where the proportionality coefficient B(R; E) depends on 
the details of the coupling. Thus finally we find 

1 - \S\\{E)\ 2 = 7r|B(R; E)\ 2 C{E), (5) 

where C{E) is the LDOS of the unperturbed system. To 
take into account the perturbation by the measurement 
antenna, in the correct derivation the Green's function 
G has to be replaced by a properly renormalized Green's 
function £p, where /3 corresponds to the scattering length 
of the antennai 18 ' 19 The coefficient |-B(R; E)\ 2 is denned 
by the antenna properties which are slowly varying with 
energy and thus can be set to \B(R; Eo)\ 2 in the vicinity 
of the Dirac point E D = k 2 D . With C(E) ~ \E - E D \ ~ 
\k — ku\ we expect 

1 - \S U (E)\ 2 ~ \k - k D \. (6) 

Now we can proceed to the discussion of our exper- 
imental findings for the graphene layers. We generated 
different graphene flakes and averaged over flakes and an- 
tenna positions. By this averaging procedure the LDOS 
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FIG. 5. Spectrum of a graphene flake, where the antenna is 
placed at the zigzag boundary. The inset shows the actual 
lattice and the cross marks the position of the measuring an- 
tenna. The dashed vertical line is the resonance position of a 
single disc. 

is turned into the DOS. The result is shown in Fig. 2) 
At the frequency corresponding to the eigenfrequency a 
single disk (vertical dotted line) a spectral gap is ob- 
served which corresponds to the Dirac point in an infi- 
nite system. Additionally a symmetric linear increase in 
the DOS is found close to the Dirac point. The overall 
decrease in 1 — |Sii| 2 with frequency has its origin in the 
fact that the measuring antenna is not placed at a disk 
center, but somewhere in between the discs, leading to 
a frequency dependent decrease in B at the antenna due 
to Eq. ©. 

Now let us have a look what happens if we move the 
antenna to the outside, e.g. to the zigzag edge. This 
measurement is presented in Fig. [S] A state close to the 
Dirac point is observed, which cannot be seen in Fig. |4l 
even though it is the same system. As the measured 
quantity is related to the LDOS and not the DOS it is 
possible to test the states locally and not only globally. 
We also moved the antenna to the armchair edge and 
did not observe a state close to the Dirac point. This 
is in accordance with the findings for a zigzag-edged in- 
finitely long graphene ribbons where the effect of edge 
states was discussed^^ and where the localization of the 
edge states was estimated to about one to two atomic lay- 
ers. If we move the antenna to the lower left corner (see 
inset Fig. [5]) where there is an isolated disk not contained 
in a full benzene ring, again a single state is observed at 
the Dirac point. By removing the disk at the corner the 
state disappears. 

To interpret our experimental findings we additionally 
performed numerical calculations by determining eigen- 
values and eigenvectors of the Hamiltonian 

H = E D • 1 + kC, (7) 

where Ed is the gap energy, k the coupling constant, and 
C the connectivity matrix with elements C nm — 1 or 
for nearest and non-nearest neighbors, respectively. The 
resulting DOS is shown in Fig. [BJ The histogram corre- 
sponds to the experimental findings presented in Fig. (SJ 
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FIG. 6. (color online) DOS from a numerical calculation of 
the systems shown in the top right corner (histogram). Ad- 
ditionally the right inset shows an grey scale indication of an 
edge state found close to the Dirac frequency. In the left inset 
an typical eigenfunction of an edge state is presented, where 
the disk in the lower left corner was removed. The crosses cor- 
respond to a system with periodic boundary conditions taking 
178 disks into account. 

where the edge states are strongly observed but inter- 
nal states are suppressed. One observes a peak of the 
DOS at the Dirac point. Analyzing the wave functions 
at this energy we found only edge states similar to the 
one presented in the left inset of Fig. [6] as well as a single 
corner state observed if an individual disk not part of a 
full benzene cell is added (see right inset). By moving the 
measurement antenna toward the center of the system, 
corner or edge states cannot be seen as their probability 
inside is very small and the measured quantity is related 
to the LDOS. In addition we performed the calculation 
with periodic boundary conditions. Now the edge states 
have disappeared, and the triangular gap at the Dirac 
point is clearly seen (crosses in Fig. [5J). This is similar to 
the situation of an internal antenna, i.e., it corresponds 
to Fig. |4l where the Dirac point is found as well. In the 
numerics with periodic boundary conditions edge states 



can not exist, whereas in the experiment the existing edge 
states are not observable, due to their small amplitude 
inside the flake. Finally by removing the corner disk 
leaving the lattice only with complete benzene cells the 
corner state did vanished, again in accordance with the 
experimental findings. 

In this paper we presented the first microwave real- 
ization of graphene in the tight-binding approximation. 
Using high index of refraction disks inside two metal- 
lic plates, where the disks resonances are only coupled 
evanescently, a tight-binding realization became possi- 
ble. The Dirac point and a linear increase in the DOS 
close to it have been found as well as edge and corner 
states. It was essential for this approach to conceive a 
system where each disk brings in just one state coupling 
evanescently with its neighbors. This is a qualitatively 
new feature for the microwave experiments as a whole, 
enabling a new class of experiments in the context of 
molecular orbital theory, going by far beyond the specific 
example of graphene studied in the present work. The 
results for benzene shown in Fig. [3] illustrate for a text- 
book example the feasibility of the method. Due to the 
flexibility of the set-up and the closeness of the system 
additionally graphene like systems with disorder, defects, 
graphene quantum dots^i or even Dirac oscillators^ can 
be investigated. By introducing different types of disks 
with different resonance frequencies it is also possible to 
investigate systems with different atom per unit celllike 
boron-nitride, corresponding to a Dirac equation for par- 
ticles with masses. 
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